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The resonance parameters of the excited w-family resonances, namely the w(4040), (4160), and 
(4415), were determined by fitting the R-values measured by experiments. It is found that the 
previously reported leptonic partial widths of these states were merely one possible solutions among 
a four-fold ambiguity. By fitting the most precise experimental data on the R-values measured by 
the BES collaboration, this work presents all four sets of solutions. These results may affect the 
interpretation of the charmonium and charmonium-like states above 4 GeV/c’. 


PACS numbers: 12.38.Qk, 13.66.Bc, 14.40.Gx 
I. INTRODUCTION 


Charmonium, bound state of a charm and an anti- 
charm quarks, is one of the most interesting two-body 
systems which were studied extensively in particle physic- 
s. Although the first charmonium state was discovered 
more than thirty-five years ago, there are still many 
puzzles in charmonium physics. The charmonium spec- 
troscopy below the open charm threshold has been well 
measured and agrees with the theoretical expectation- 
s (such as potential models and lattice QCD); howev- 
er, for the charmonium states above the open charm 
threshold, there are still lack of adequate experimental 
information and solid theoretical inductions. For exam- 
ple, recently many new particles have been discovered, 
named XY Z-particles, and the overwhelming vector s- 
tates in the 4 GeV/c? to 5 GeV/c? mass range make the 
classification of these states as the charmonia question- 
able [1-3]. In explaining these vector charmonium states, 
the leptonic partial widths provide very important infor- 
mation. As we know, the vector quarkonium states could 
be either S-wave or D-wave spin-triplet states, with the 
S-wave states couple strongly to lepton pair while the D- 
wave states couple weakly since the latter are only pro- 
portional to the second derivative of the wave-function 
at the origin squared, as expected in the potential mod- 
els. This leads people believe that the (4040) is the 3S 
charmonium state, w(4160) the 2D state, and (4415) 
the 4S state. This has been a well accepted picture 
for more than two decades before the discovery of the 
so-called Y particles, namely, the Y(4008) and Y (4260) 
observed in ete~ — ata J/y final state [4, 5], and 
the Y (4360) and Y (4660) observed in ete~ > ntry’ 
final state [6, 7]. With seven states observed between 
4.0 GeV/c? and 4.7 GeV/c?, some people started to cat- 
egorize some of these as non-conventional quarkonium 
states, while others tried to accommodate all of them in 
modified potential models. Many of the theoretical mod- 
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els use the leptonic partial widths of these states to dis- 
tinguish them between S- and D-wave assumptions [8, 9], 
and most of the time, the values on the leptonic partial 
widths are cited from the PDG [10] directly. Although 
the resonance parameters of these excited w states have 
been measured by many experimental groups, all of them 
were obtained by fitting the R-values measured in the rel- 
evant energy region. The most recent ones, which were 
from a sophisticated fit to the most precise R-values mea- 
sured by the BES collaboration [11, 12], are the only 
source of the leptonic partial widths of these three w s- 
tates now quoted by the PDG [10]. 

In fitting to the BES data, unlike the previous anal- 
yses, the BES collaboration considered the interference 
between the three resonances decaying into the same fi- 
nal modes, and introduced a free relative phase for the 
amplitude of each resonance [13]. The new parametriza- 
tion of the hadronic cross section results in a prononce 
increase of the (4160) mass, and significant decrease of 
the leptonic partial widths of (4160) and (4415). 

As has been pointed out in a recent study [14], there 
are multiple solutions in fitting one dimensional distribu- 
tion with the coherent sum of several amplitudes and free 
relative phase between them. Exactly the same kind of fit 
to the R-value is performed in the present study, multi- 
ple solutions are indeed found in extracting the resonance 
parameters of the excited w states. The only difference 
between these multiple solutions is the coupling to the 
e+e, namely the leptonic partial width, while the mass- 
es and the widths of the resonances remain the same for 
all the solutions. 

In the following, we firstly introduce a simplified fit 
scheme similar to that used by the BES collaboration, 
and extract the resonance parameters. Then we study 
the multiple solution problem in the light of the toy sim- 
ulated data for illustrative purpose. We will discuss the 
consequence of the multiple solutions in fitting the R- 
values at the end of the paper. 


II. THE FIT TO THE R DATA 


To faciliate our study, only the R-values provided by 
the BES collaboration are used [11, 12], and in the data 


fitting, only statistical uncertainties are considered, as 
the systematic errors at all the energy points are highly 
correlated. 

We fit the ete” annihilation cross section o(ete~ ) = 
R - 86.85/s (s in GeV? and ø in nb). The standard 
chisquare estimator is constructed as 
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where o; ™ and Ao;*”’ indicate respectively the experi- 
mentally measured cross section and its error at the l-th 
energy point (the number of points is denoted as Npt), 
while gihe is the corresponding theoretical expectation, 
which is composed of two parts 


ate (s) = a”®® (s) 4 geo (s) 3 (2) 


Here o°°”": is the contribution from continuum and is pa- 


rameterized simply as 


a”: (s) — A+ Bis — 2Mp+) ; (3) 


where A and B are free parameters, and Mp+ is the mass 
of the charged D meson. g”®°®: is the contribution from 
the resonances. Here following previous analyses [13, 15], 
the assumption that the continuum production and the 
resonance decays don’t interfere with each other is adopt- 
ed. The three wide resonances shown in the data are 
close and have the same decay modes, the interferences 


between them must be included. Then the amplitude 


reads 
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with T$, r}, rs°, and M; denoting total width, partial 
width to hadrons, partial width to et e~ pair, and mass of 
the resonance j, respectively. The total amplitude, which 
is the coherent sum of the three resonances, once squared, 
contains the interferences of the type RTT}. Here ġ; is 
a phase associated to resonance j. If the resonances are 
quite broad, the interference effect will distort the shape 
of the resonances, the width might appear broader or 
narrower, and the position of the peak might be shifted 
as well. 
The total cross section of the resonances 
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Since what we actually obtain is the squared module of 
the amplitude, only two relative phases are relevant. 

By minimizing the y? defined in Eq.(1), we obtain the 
results as displayed in Fig. 1. Just as expected, the inter- 
ference effect changes the shape of each resonance signifi- 
cantly: the ~(4040) becomes narrower while the ~(4160) 
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FIG. 1: Total hadronic cross section in nb obtained as 


o(ete~ — hadrons) = R- 86.85/s (s in GeV”) from R mea- 
surements of BES [11, 12] and the fit to data. The curves 
show the best fit (identical for all the four solutions), the con- 
tribution from each resonance, as well as the interference term 
(one of the four solutions). 


and y(4415) become wider than the previous published 
results where the interferences between resonances are 
neglected [15]. 

The fit results are presented in Table I, there are four 
solutions found in the fit. It should be noted that the 
four solutions have identical x7, masses, and total width- 
s for the resonances, but different partial widths to lepton 
pairs. One may suspect that the existence of four solu- 
tions is due to low precision of the measurements. We 
will show in the next session that multiple solution is a 
real effect. The improvement of the precision of the mea- 
surements can not change the fact that there are four 
solutions in fitting these data. 


II. MULTIPLE SOLUTIONS 


The existence of four solutions in the fit described 
above can be tested through the toy simulated data. The 
special steps for this approach is: 


1. The curve of the best fit in previous session (with 
background contribution subtracted) is used as the 
probability density function (PDF) to compute 
a set of experimental points at different energies 
(Np: = 100 in this work) within the range 3.75 GeV 
< Js < 4.85 GeV. 


TABLE I: Four groups of solutions for the data fitting. The 
four solutions have identical resonance masses (M) and total 
widths (T+), but significant different leptonic partial widths 
(Tee) and the relative phases (¢). The fit yields x? = 91, 
A = (15.05 + 0.59) nb, and B = (—1.64 + 0.67) nb/GeV for 
all the solutions. 


Parameter (4040) (4160) (4415) 
M (MeV) 403446 419347 4412415 
T; (MeV) 87+11 79+14 118432 
TY) (keV) 0.66+0.22 0.42+0.16 0.45+0.13 
o (radian) 0 (fixed 2.7+0.8 2.00.9 
TS (keV) 0.72+0.24 0.7340.18 0.60+0.25 
¢ (radian) 0 (fixed 3140.7 1.41.2 
TS? (keV) 1.28+0.45 0.62+£0.30 0.59 + 0.20 
¢ (radian) 0 (fixed 3.7+0.4 3.8+0.8 
TO (keV) 1.414012 1.1040.15 0.78 0.17 
o™ (radian) 0 (fixed 41+0.1 3.2403 


TABLE II: Fit results for four groups of solutions with the 
toy simulated data points. The definitions of the parameters 
are the same as in Table I. The fit yields x? = 1.0 x 107° for 
all the solutions with the expectation of y? = 0. 


Parameter (4040) (4160) (4415) 

M (MeV) 4033.5+0.3 4192.8+0.3 44124+04 
r+ (MeV)  87.2340.49 79.0040.53 118.11 + 0.56 
TP (keV) 0.6644£0.005 0.417+0.004 0.454 + 0.003 
o (radian) O (fixed) 2.701 + 0.012 2.002 + 0.012 
T? (keV) 0.723+0.006 0.73140.005 0.596 + 0.003 
® (radian) 0 (fixed) 3.051 0.001 1.432 + 0.014 
TS) (keV) 1.283 0.005 0.620 0.006 0.590 + 0.003 
® (radian) 0 (fixed) 3.732 0.006 3.789 0.013 
TÉ (keV) 1.397 0.006 1.087 0.008 0.774 0.003 
p% (radian) O (fixed) 4.082 0.005 3.218 + 0.009 


2. A relative error of 1% and an absolute error of 
0.01 nb are added in quadrature and the total error 
is assigned to each data point. The inclusion of a 
0.01 nb absolute error is to weaken the chisqure- 
weight of points with small cross sections. 


Fit the simulated data with the similar y? in E- 
q. (1) but with only resonance cross section included (i.e., 
othe-(s) = o"®® (s) in Eq. (2)), the four groups of solu- 
tions obtained are summarized in Table II and displayed 
in Fig. 2. Comparison of two tables obviously indicates 
that the central values of the parameters are consistent 
with each other. 

From Table II and Fig. 2, we can see that the largest 
Pee is more than twice the smallest value in the four so- 
lutions. We also notice that the first solution is identical 
to the BES published one [13], and the second solution is 
about the same as the one listed as the best estimation 
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FIG. 2: Four groups of solutions obtained from the fit to the 
generated data. For clearness, the interference curves have 
been moved downward by 10 nb, the black dashed line at 
—10 nb corresponds to zero cross section. 


of the partial widths of these states by the PDG [10]. 


IV. DISCUSSIONS 


As the Pee of the vector resonances are closely related 
to the nature of these states, the choice among the dis- 
tinctive solutions affects the classification of the charmo- 
nium and charmonium-like states observed in this energy 
region. 

The calculation of the ete~ partial widths of the S- 
wave charmonium states is well summarized recently in 
Table 2 of Ref. [16]. We can see clearly that many of 
the theoretical calculations give large w(3S) and 7(4S) 
decay widths compared to the PDG values (about the 
same as the second solution listed in Table I above), but 
the agreement with the third or the fourth solution is 
much better. 

The Y (4260) was proposed to be the w(4S') state 
and the (4415) be ~(5S') in Refs. [8, 9], we can see 
that in this assignment, the calculated partial widths of 
w(4040) = ¥(3S) and (4415) = (5S) [9] agree well 
with the fourth solution listed in Table I. 

Of course the possible mixing between S- and D-wave 
states will change significantly the theoretical predic- 
tions of the partial widths of these states [17], and the 
QCD correction, which is not well handled [16], may also 
change the theoretical predictions significantly. So far, 
we have no concrete criteria to choose any one of the 


solutions as the physics one. 

It should be noticed that if the Y states are considered 
together with the excited ~ states in fitting the R-values, 
there could be even more solutions, the situation may 
become more complicated. 

We also notice that the existence of the multiple solu- 
tions is due to the inclusion of a free phase between two 
resonances, if these phases can be determined by other 
means (either theoretically or experimentally), it will be 
very helpful to know which solution corresponds to the 
real physics. 


V. SUMMARY 


Based on the R scan data, the resonance parameters 
of excited w-family resonances are fitted. We found that 
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